KREIN'S RESOLVENT FORMULA FOR 
SELF-ADJOINT EXTENSIONS OF SYMMETRIC 
SECOND ORDER ELLIPTIC DIFFERENTIAL 
OPERATORS 

ANDREA POSILICANO AND LUCA RAIMONDI 

Abstract. Given a symmetric, semi-bounded, second order ellip- 
tic differential operator A on a bounded domain with C^'^ bound- 
ary, we provide a Krem-type formula for the resolvent difference 
between its Friedrichs extension and an arbitrary self-adjoint one. 



1. Introduction. 

Given a bounded open set ^2 C M", n > 1, let us consider a second 
order elliptic differential operator 

n n 
i,j=l i=l 

Such an operator A, under appropriate hypotheses on its coefficients 
and on Q (these will be made precise in the section 3), is closable and 
its closure A^[^, the minimal realization of A, has domain given by 
iJ^(n), the closure of C^{n) with respect to H^{n) Sobolev norm. If 
A is symmetric then Amin is symmetric but not self-adjoint, i.e. A is not 
essentially self-adjoint. Indeed ^min = ^max, where Amax, the maximal 
realization of A, has domain made by the functions u G L^(f2) such 
that Au G L^(f2). Assuming that A^i^ is semibounded, then Amin has 
a self-adjoint extension Aq (the Friedrichs extension, corresponding to 
Dirichlet boundary conditions), Amin Q Aq C Amax, and hence Amin 
has infinitely many self-adjoint extensions. 

The problem of the parametrization of all self-adjoint extensions of 
Amin in terms of boundary conditions was completely solved (in the case 
of an elliptic differential operator of arbitrary order) in [12] (for some 
older papers about similar topics we just quote [5] and [23]). Here, 
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by using the approach developed in [T6]-[l9j, we give an alternative 
derivation of such a result by providing a Krem-like formula for the 
resolvent difference between an arbitrary self-adjoint extension of Amin 
and its Friedrichs extension Aq. For the sake of simplicity here we 
consider the case of a second order differential operator. The case of 
higher order operators can be treated in a similar way. 

In the case A is the Laplacian, the Krem resolvent formula here 
presented has been given in [TH], Example 5.5. For other recent results 
on Krem- type formula for partial differential operators see [I], ^22j, [4j, 
i, i, [3]. 

In order to help the reader's intuition on the results here presented, 
in Section 4 we consider one of the simplest possible examples: a 
rotation-invariant elliptic operators A on the disc D C M^. Thus, 
notwithstanding the symmetric operator here considered has infinite 
deficiency indices, due to the presence of symmetries the resolvents of 
their self-adjoint extensions can be written, by separation of variables, 
in a form which resembles the finite indices case (see the comments in 
Remark 4.1), and the corresponding spectral analysis becomes simpler. 
As illustration, given any sequence {A„}^ C M, boundary conditions 
at dD can be given for which A is self-adjoint and such that {A„}^ is 
contained in its point spectrum. Remark shows that such boundary 
conditions can be quite different from the usual ones. 

2. Preliminaires 

For the reader's convenience in this section we collect some results 
from [TB]-[T1]. We refer to these papers, in particular to [T^], for a 
through discussion about the connection of the approach here presented 
with both the standard von Neumann's theory of self-adjoint extension 
[15] and with Boundary Triple Theory [6]. [TO]. 

From now on we will denote by 

^(L), ^(L), p(L) 

the domain, kernel, range and resolvent set of a linear operator L. 
Let be a Hilbert space with scalar product (■, ■) and let 

Ao : ^{Ao) CJif^Ji^ 

a self-adjoint operator on it. We denote by J^Ao be the Hilbert space 
given by the linear space 3i{Aq) endowed with the scalar product 
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Given then a Hilbert space f) with scalar product (■,■) and a hnear, 
bounded and surjective operator 

such that ^{t) is dense in , we denote by S the densely defined 
closed symmetric operator 

S : ^(r) C ^ ^ , 50 := Ao0 . 

Our aim is to provide, together with their resolvents, all self-adjoint 
extension of S. 

For any z G p{Aq) we define the bounded operators 

R,:={-A^ + z)-^ -..^ ^J^A,, 

(2.1) :=(ri?,)*:f)^jr. 

By [IT], Lemma 2.1, given the surjectivity hypothesis ^{t) = f), the 
density assumption ^{t) = is equivalent to 

MiG,) n &{Ao) = {0} . 

However, since by first resolvent identity 

(2.2) (z - w)R^G, = G^-G,, 
one has 

^{G^ - G,) C ^(Ao) . 

Prom now on, even if this hypothesis can be avoided (see [I6]-[T9]), for 
the sake of simplicity we suppose that 

G p{Ao) . 

We define the family F^, z G p{Aq), of bounded linear maps 

(2.3) F, : f) 1^ , F, := r (Go - G,) = -ztA.'G, . 
Given then an orthogonal projection 

n : [) ^ [) , 1^0 = m^) , 

and a self-adjoint operator 

e : ^(6) C 1^0 t)o, 
we define the closed operator 

r.,n,e := (e + n F, n) : ^(9) C [)o ^ f)o, 
and the open set 

Zu,e ■■= {z G p(v4o) : G p(F,,n,e)}. 
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With such premises the next two theorems have straightforward proofs. 
Theorem 12.11 is an obvious modification (taking into account the hy- 
pothesis G p{Aq)) of Theorem 3.1 in [18] (also see [17], Theorem 3.4); 
Theorem 12.21 is the combination of Theorem 12. II with Theorem 2.1 and 
Theorem 2.4 in [19] (also see [16], Theorem 2.1, [17], Theorem 2.2, for 
the case 11 = 1). 

Theorem 2.1. The adjoint of S is given by 

S* : ^{S*) CJif-^J^, S*(f) = Ao0o , 

^{S*) = {0 e ^ : = 00 + Go(^, 00 e ^{Ao), (^ei)}. 
Moreover 

(2.4) V0 , ^ G ^{S*) , {S*4>, ^) - (0, S*ij) = (r0o, C^.) - iC<t>, r^Po) ■ 
Theorem 2.2. The set Zu,0 is not void, 

C \ M C Zn,e 

and 

Rz,u,e '■= Rz + c^n Q n , z e Zu^e , 

is the resolvent of the self-adjoint extension Ah q of S defined by 
Au,e : ^{An,e) C ^ , An,©0 = 5*0 = Ao0o , 
^(An,e) = {0 G ^{S*) : G ^(0) , nr0o = 6(4 . 
Remark 2.3. Note that, since 0o = ^5'*0, 

nr0o = ec^ ^ nro0 = ec<^ , 

where the regularized trace operator fo is defined by 
To : &{S*) ^ i) , fo0:=rAo'5*0. 

By exploiting the connection with von Neumann's theory (see [19], 
section 3; see also [17], section 4, for the case of relatively prime exten- 
sions) one obtains 

Theorem 2.4. The set of operators provided by Theorem \2.S\ coincides 
with the set S'{S) of all self-adjoint extensions of the symmetric oper- 
ator S. Thus S{S) is parametrised by the bundle p : E(J)) — > P(f)), 
where P(P)) denotes the set of orthogonal projections in [) and p'^(Jl) 
is the set of self-adjoint operators in the Hilbert space ^{IV). The set 
of self-adjoint operators in (), i.e. p~^{l), parametrises all relatively 
prime extensions of S i.e. those for which ^{A) n ^{Aq) = S>{S). 

We conclude this section with a result about the spectral properties 
of the extensions (see [6j, Section 2, for point 1 and ^18j, Theorem 3.4, 
for point 2). 
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Theorem 2.5. 1) 

A G (Tp{An,e) n p{Ao) G (rp(Tx,u,e), 

where crp(-) denotes point spectrum. An analogous result holds for the 

continuous spectrum. 

2) 

Gx : ^(rA,n,0) ^ ^i-Au,e + A) 
is a bijection for any A G crp(An,0) H p(Ao)- 

3. Extensions and Krein's Formula. 

Let f2 C M", n > 1, a bounded open set with a Lipschitz boundary. 
We denote by H''{fl) the Sobolev-Hilbert space given by closure of 
C°°{Cl) with respect to the norm 

0<aiH \-an<k 

Analogously -ffo(^) £ H''{il) denotes the closure of C^(fi) with re- 
spect to the same norm. 

Given the differential expression 

n n 

A = V-aV — 6-V — c = ^ di{aijdj) — ''^hdi — c 

i,j=l i=l 

we suppose that the matrix a{x) = {aij{x)) is Hermitean for a.e. x & Q, 
that there exist yUi > 0, /i2 > such that 

V^GM", f^im'<^-a{x)^<fi2Uf 

and that 

6i G , c G , q = n if n>3,q>2 a n = 2. 

Then A maps H'^^Q) into H-^{n) (see e.g. [7], Section 1, Chapter VI), 
where H~^{Q) denotes the adjoint space of Hq{Q), the sesquilinear 
form 

qA : X H^{Q) ^ C 

qA{u,v) := - ((VM,aVt;)L2(f^) + {u,b-Vv) l2(^^) + {u,cv)l2^q)) , 

is continuous and there exists a positive constant A such that — + A 
is coercive (see e.g. Proposition 1.2, Chapter VI). Thus by Lax- 
Milgram Theorem (see e.g. [7], Theorem 1.4, Chapter VI) there exists 
a unique closed, densely defined, linear operator 

Ao : ^{Ao) C L^in) L^in) , A^u = Au , 

&{Ao) = {ue Hl{n) : Au G L\Vl)} , 
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such that 

Vm G ^(Ao) , \/v G H^in) , QAiu, v) = {u, Aov)L2(n) . 

Moreover S>{Aq) is dense in HI{VL), G p{—Aq + A), Aq has a compact 
resolvent and its spectrum consists of an infinite sequence of eigenvalues 
each having finite multiplicity and with ReA„ < —A. An analogous 
result holds for the sesquilinear form 

qXiu.v) := qA{v,u) 

and the operator corresponding to q\ is the adjoint A^. 
Suppose now that 

diQij G L'^{Vl) , g = ?T,if?T,>3,g>2ifn = 2, 

so that, by Sobolev Embedding Theorem, A is continuous from H'^iVt) 
into L^(fi) and 

By interior regularity estimates (see e.g. ||13|, Section 7, Chapter 3) 
ylfc'^(n), the restriction of A to C^(f2), is closable and its closure is 
given by Amin C Aq, the minimal realization of A, defined by 

Prom now on we suppose that 

(1A = (1*A- 

Thus Aq is a self-adjoint operator, the Friedrichs extension of the closed 
symmetric operator Amin and one has 

^min = {Mc^{n)Y = Amax , 

where Amax, the maximal realization of A, is defined by 

max^ ■ — Au , 

^(Amax) := G : Au G ^^(l])} . 

Hence 

^(Ao)=i/o'(^)n^(Amax). 

Moreover 

SO that A\c^{n) is not essentially self-adjoint, 

A . C Ar. C A 

^min ^ ^0 -5: ^max ; 

and the symmetric operator Ajam has infinitely many self-adjoint ex- 
tensions. We want now to find all such extensions and to give their 
resolvents. In order to render straightforward the application of the 
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results given in Section 2, we would like to have a more explicit char- 
acterization of !^{Aq). Thus in the following we impose more stringent 
hypotheses on the set fl. 

Suppose that the boundary of f2 is a piecewise surface with cur- 
vature bounded from above and that atj G C{Q) when n > 3. Then, 
by global regularity results (see e.g. [I3], Chapter 3, Section 11), the 
graph norm of Amax is equivalent to that of H'^{fl) on C^{Cl), the space 
of smooth functions on Q which vanish on its boundary dQ. Thus 
^\cg°{n)y the restriction of A to C^{Q), is closable and its closure is 
given by 

io : ^0 (fi) C L^{n) L^{VL) , Aqu := Au , 

where Hq{Q) denotes the closure of C^{Cl) with respect to the H'^{Q) 
norm. 

Without further hypotheses on Q, Aq ^ Aq is possible: for example 
if n is a non-convex plane polygon then the Laplace operator A is not 
self-adjoint on Hq{Q). Indeed by [2] it has deficiency indices {d-, d+) = 
{d, d), where d is the number of non-convex corners . 

Suppose now that the Lipschitz continuous up to the bound- 

ary and that dQ is C^'^, i.e. it is locally the graph of a function 
with Lipschitz derivatives (see e.g |Tl], Section 1.2, for the precise def- 
inition). Then (see e.g. [H], Chapter 1, Section 8.2, [TI], Section 1.5) 
there are unique continuous and surjective linear maps 

p : H\n) H^/\d^) , 

7„ : H\n) ^ H'^'^dn) ® H'/\dQ) , 7,0 := (p0, r,0) , 
such that 

pcf) (x) := (x) , Ta(t> (x) = — (x) := ^ aij{x)ui{x)dj(j) (x) 

for any G C°°{Q) and x G dQ. Here u = {ui, . . . , i/„) denotes the out- 
ward normal vector on dQ and H^{dQ), s > 0, are the usual fractional 
Sobolev-Hilbert spaces on dfl (see e.g. [H], Section 1.3.3). Moreover 
Green's formula holds: for any u G H'^{Q) and v G H^{VL) fl HI{VL) one 
has 

(3.1) {Au, v)l^(q) = {u, Aov)L2^n) - {pu, ^a^^)L2(an) ■ 

By proceeding as in the proof of Theorem 6.5 in [H], Chapter 6 (which 
uses (13. ip ) the map 7, can be extended to (see [TT], Theorem 1.5.3.4) 
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where H~'^{dQ) denotes the adjoint space of H^{dQ), and Green's for- 
mula fl3.ip can be extended to the case in which u G ^{A^a.^): 

(3.2) (AmaxM, v)L2{n) = {U, AoV)L2{n) - {pU, TaV)_i^i . 

Here i denotes the duahty between H^/'^{dVt) and H~~^/'^{dVt). 

With such definitions of p and r one has (see e.g. [11], Corollary 
1.5.1.6), 

Moreover, by the stated properties of p and p, by the equivalence of 
the graph norm of Amax with the H'^{Q) norm on Hq{Q) and by the 
density of C°°{Cl) in ^(Amax), one gets the equalities 

^o'(^^) = H^{^) n H^{n) = ^(A^ax) n H^{n) = ^{Ao) , 

so that Aq = Aq. 

In conclusion we can apply the results given in Section 2 (by adding, 
if necessary, a constant to Aq we may suppose that G p{Aq)) to the 
self-adjoint operator 

Aq : H^{n) n H^{n) c L'^{n) , Aqu ■.= Au , 

with S = A^in, f) = H^/\dn) and 

r : H\Q) n //o^(fi) ^ H^/\dQ) , r := r„ 

Note that jr(r) = i/o(^) since ^(7^) = ^(71) by z/(x)-a(x)z/(x) > 
/xi > 0, X G c)f2, and that r is surjective by the surjectivity of 7a. 

Thus, by Theorem 12. 4[ under the hypotheses above, the set i^(Amin) 
of all self-adjoint extensions of A-^^^ can be parametrized by the bundle 

Now, in order to write down the extensions of A^[^ together with their 
resolvents, we make explicit the operator Gz defined in f l2.ll) . By The- 
orem [5TT], since Amax = ^min) havc 

^(Amax) = {U = U0 + G^h , Uq G H\^1) f] H^{n) , h G H^/^{dn)} , 

^maxM = ^0^0 • 

Thus A^^^Goh = and so by there follows, for all h G H^''^{d9) 
and for all u G ^{Aq), 

{Goh, AqU) VHn) = {pGoh,TaU)^i l . 

Since, by ([MD, 

{Goh,AQu)L'2(n) = (Co/i, v4inaxM)L2(Q) 
= {Goh,A*^^^u)L'2(n) = {h,TaU)Hi/2^gQ) , 
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one obtains pG^h = Ah, where 

A : H^/^{dn) ^ H-^/^{dn) 

is the unitary operator defined by 

For successive notational convenience we pose S := A~^. 

Remark 3.1. If dQ carries a Riemannian structure then H'^{dfl) can 
be defined as the completion of C°°{dQ) with respect of the scalar 
product 

{f,g)H'>{dn) ■■= if, {-^LB + iyg)L^{dn) ■ 

Here the self-adjoint operator A^^ is the Laplace-Beltrami operator in 
Li^^dVt). With such a definition (—Alb + 1)^^^ can be extended to the 
unitary map A. 

Since Gz = Gq + zA^^Gz by (12. 2p . Gzh is the solution of the Dirichlet 
boundary value problem 

„N J ^maxG^z^ = zGzh , 

^ ■ ^ \pGzh = Ah. 

Thus we can write GqS = isT, where K : H-^/'^{dn) ^{A^^^) 
is the Poisson operator with provides the solution of the Dirichlet 
problem with boundary data in H~^/'^{dVL). Analogously we define 
Kz : H~^/\dn) ^(A^ax) by Kz := G^S. Note that Goh, hence 
Gzh, is uniquely defined as the solution of (13. 3p : for any other solution 
u one has u — G^h G .J^{Aq) = {0}. 

Now, according to (12. 3p . we define the bounded linear operator 

r, : H'/'^idn) ^ H'/\dn) , r, := r(Go - Gz) , 

which, by (12. 2p and the definitions of K and Kz, can be re- written as 

(3.4) = -ZTaA^'Gz = Z TaRzKA = (faK - f„K,)A . 

By pGoh = Ah, by Theorem 12.11 and by Remark 12.31 we can define the 
regularized trace operator 

fafi : ^(Amax) ^ H^^^dn) , 

(3.5) fa^oM := Taiu - GoT,pu) = faU - PaPU = TaA^^ A^^y^U , 

where the linear operator Pa, known as the Dirichlet-to- Neumann op- 
erator over dQ, is defined by 

Pa : H-'/\dn) ^ H-'/''idn) , Pa:=faK. 

In conclusion, by Theorems 12.21 and 12.41 one has the following 
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Theorem 3.2. Any self-adjoint extension A o/Amin is of the kind 

A : ^(i) c L\n) L\n) , Au = a^^^u , 

^(i) = {ue ^(A^ax) : Spn e ^(9) , Ufa,oU = eSpu} , 

where {U,Q) e E{H^/'^{dn)), and 

i-A + z)-^ = i-Ao + z)-' + G,U (6 + n T,U)-^UGl , 

with Tafl, Gz and defined by hS. 5\) . h3. 3\) and \3.4\) respectively. 

Remark 3.3. By preceding as in [19j, Example 5.5, in the case the 
L^(9f2)-symmetric, bounded linear operator _B : H^/'^{dVL) — > H^^'^{d^l) 
is such that 6^ := (-P„ + 5)A, ^{Qb) = H^^^{dn), is self-adjoint {B 
pseudo-differential of order strictly less than one suffices), the exten- 
sion Ab corresponding to (1,9^) has domain defined by Robin- type 
boundary conditions: 

^{Ab) := {u G H\n) : r^u = Bp} . 
4. A SIMPLE EXAMPLE. 

One of the simplest examples is given by a rotation invariant 2nd 
order elliptic differential operator on the unit disc D C M^. Thus we 
consider the self-adjoint extensions of 

A^in : H^{D) C L\D) ^ L\D) , A^^^u = Au 

where 

y4 = V-aV — c, aij{x) = a {\\x\\) 6ij , c{x) = c {\\x\\) . 

We suppose that a is Lipschitz continuous, info<r<i ci{r) > 0, and that 
c G L'^{{0, l);rdr), q > 2. By adding, if necessary, a constant to c we 
suppose that —Aq > 0. 

In L'^{D) ~ -^^^((0, 1); rdr) (g) L'^{{0, 27r); d(p) we use the orthonormal 
basis {Umn}, m E N, n E Z, 




made by the normalized eigenf unctions of the Friedrichs extension Aq of 
A. Here {Umn}, m eN, is the orthonormal basis in -Z^^((0, l);rdr) made 
by the normalized eigenfunctions of the self-adjoint Sturm-Liouville 
operator 

Lnfir) = -- irair)f'ir))' + (c(r) + ^"j /(r) , n > , 
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with boundary conditions = and limr^o+ ^f'{f) = if n = 0, 

/(0+) = if n 7^ 0. Denoting by > 0, m E N, the eigenvalues of 
Ln, one has 

a{Ao) = (TdiAo) = {-A^|„| , meN,neZ}. 
In 7/^/^(5^) we use the orthonormal basis {ck}, /c G Z, defined by 

V^(P + 1)1/4 • 

We want now to compute the matrix elements, relative to the basis 
{Umn}, of the resolvents of the self-adjoint extensions of A^^i^. 
By defining 

u^n ■■= lim a(r) M^„(r) , 

r Jl 

one has 



[Cojmnfe '■ — {Umn, GoCk) L^(D) — {Gyjmni ^k) H^/'i{S'^) —'■ [G}^kmn 

^m\n\ 



m\n 

Since = Gq — z{—Ao + z)'^Go, one has then 

[G^jmnfc [G^J/jmn [GojrriTifc T2 ~j~ [Gojmnfc 

■^m\n\ + ^ 



m\n\ m\n\ 

Analogously 

[r^jifc := - z{ei, Ta{-AQ + ^)~"^Goefc)j:^i/2(5i) 

00 2 

^-.(P + i)'/^^ a""" 

m=l '^mlfcll^mlfcl + 

Thus, in the case the orthogonal projection 11 is the one corresponding 
to the subspace of H^^'^{S^) generated by {e^ , k E I}, I C Z, and 
[0]ifc = ^k^ik, A; G /, by Theorem 12.21 one obtains 

(-^n,e + z)'-^ := {Umn, (-^n,e + z)~'^Ufnn) l'^(d) 

J mnfhh 



for any n E I, and 



-An,e + z)-^ 



■'mm'-'nn 
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for any n ^ I. Once the resolvent has been written as above, by 
Theorem 12.51 given any sequence 



posing 
one obtains 
Moreover 



{A„}„g/ C Mnp(Ao), 
C ap{An,e) ■ 



[/ = V u 



is eigenf unction with eigenvalue A„. 

Remark 4.1. The previous example can be re-phrased in the language 
of decomposable operators (see e.g. [21], section XIII. 16): the opera- 
tor is decomposable with fibers Ao(n) = and the decompos- 
able self-adjoint extensions of A^^xa have decomposable resolvents with 
fibers given by the resolvents of the self-adjoint extensions of the fibers 
ylinin(n), which are symmetric operators with deficiency indices (1, 1). 
However this approach gives a less (than the one provided by Theorem 
13.21) explicit espression for the self-adjointness domain. 



Remark 4.2. In the case a = 1, c = 0, one has 

where J„ denotes the n-th order Bessel function, /i^n is its m-th positive 
zero and Cmn is the normalization constant. Thus 

The following remark shows that the boundary conditions corre- 
sponding to couples (n, 6) of the kind above can be quite different 
from the usual ones. 

Remark 4.3. Suppose in the previous example we take a = 1, c = 0, 
i.e. A = A, and / = {0}, Aq = 0. Then 

n : H''\s') ^ c , uf = ^r f{^) 

V 27r 7o 

and O : C ^ C is the multiplication by zero, since Fq = 0. Thus 
^(An,o) = |m G ^(Amax) : Spw = const , fi^u ((^) rf^ = o| . 
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Since A = maps constants into constants, 

{u E ^(Amax) : pu = const} = {u E H^{D) : pu = const} 
by elliptic regularity, and 

fi^ou (ifi) d(p = / [tiAq^A^s.xu]{(p) dip = / Tiu{(p)dip , 
Jo Jo 

in conclusion one has 

f Z"^'" du ^ 

^(An,o) = |m G H^{D) : pu = const , J p — ((^) = j . 
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